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Abstract 

We present the results of our analysis on charmonia {J/ip and ip') hadroproduction 
taking into account higher-order QCD effects induced by initial-state radiation in a Monte 
Carlo framework, with the colour-octet mechanism implemented in the event generation. 
We find that those colour-octet matrix elements extracted so far from Fermilab Tevatron 
data for both J/ip and ip' production have to be lowered significantly. We finally make 
predictions for charmonia production at the LHC, presenting a simple code for a fast 
simulation with PYTHIA based on the colour-octet model. 
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1 Introduction 


Likely, the Large Hadron Collider or LHC will become operative by the beginning of 
the next century, offering a wide programme of exciting possibilities in Particle Physics, 
among which the origin of mass at the electroweak scale and the possible discovery of 
supersymmetric particles, the investigation of CP violation in B mesons and detailed 
studies of top quark physics. Other interesting topics related to charm and beauty 
flavours will be covered as well benefitting of a foreseen huge statistics to be collected 
with the machine running even at “low” luminosity (~ 10 33 cm 2 s” 1 ) Jl| [0]- Following 
this line of research, we shall focus in this paper on inclusive hadroproduction of heavy 
resonances, in particular charmonia J/-0 and -0' states, whose relevance in the study of 
perturbative and non-perturbative QCD will be briefly reviewed. 

To first approximation, heavy-flavour hadroproduction can be treated in the frame¬ 
work of perturbative QCD due to the relatively large quark masses 0 and consequently 
the partonic production cross sections can be convoluted with parton distribution func¬ 
tions (PDF’s) of the colliding protons. This scheme can indeed be viewed as well 
founded since the underlying parton interaction is hard, therefore providing a rea¬ 
sonable justification for such a factorization of the production process. Higher-order 
corrections to the matrix elements of the short-distance processes may be incorporated 
likely improving the accuracy of the predicted cross sections. Even complex effects 
such as multiple gluon emission can be treated analytically according to resummation 
techniques 

Nevertheless, such an essentially analytic approach may eventually oversight the real 
complexity of the full hadronic collision and a certain amount of modifications should 
be added if a more realistic description is required. For example, parton evolution either 
at the initial or final state together with fragmentation and recombination of the beam 
remnants could play an important role leading to substantial modifications in the final 
state, even for inclusive processes. 

On the other hand, event generators based on a Monte Carlo framework represent 
an alternative philosophy of coping with the complexity of very high energy hadronic 
reactions by dividing the full problem in a set of components to be dealt separately in a 
simulation chain. The hard partonic interaction, constituting the underlying or skeleton 
process of the simulation, is dressed up with soft processes leading to hundreds of final 
particles. Certainly, such an approach is plagued with various sources of uncertainty, 
especially related to the long-distance aspects of the strong interaction dynamics, as¬ 
sociated to the assumptions and modelling of the generation. However, the impressive 
amount of collected experimental data used in a heuristic way to tune those parameters 
of the event generator, gives reliability to this method to a considerable degree provid¬ 
ing a powerful tool to interpret current data and make sensible predictions for future 
experiments. 

In this work we re-examine prompt production of heavy resonances at high-energy 
hadron colliders by means of the event generator PYTHIA ||. Undoubtedly, there is a 
widespread recognition that the study of heavy quarkonia has been playing an important 
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role over the past decades in the development of the theory of the strong interaction 
and such beneficial influence seems far from ending. Moreover, a precise understanding 
of prompt charmonia production offer several advantages from both experimental and 
theoretical points of view at the LHC. 

The fact that both J/-0 and -0' resonances can decay into a oppositely charged 
lepton pair with a respectable branching fraction (especially the former), permits its 
discrimination in a huge hadronic background. Actually, charmonia production can 
be viewed as interesting in its own right for precise tests of both perturbative and 
non-perturbative QCD, search for exotic states or as a probe of the parton content 
of the proton. Furthermore, completely different processes as for example the decay 
a golden channel to investigate CP violation in B decays, require a precise 
estimate of sources of background coming from random combinations involving prompt 
J/0’s. 

On the other hand, the experimental discovery at Fermilab 0 of an excess of in¬ 
clusive production of prompt heavy quarkonia (mainly for J/0 and 0' resonances) in 
antiproton-proton collisions triggered an intense theoretical activity beyond what was 
considered conventional wisdom until recently 0 ||. The discrepancy between the so- 
called colour-singlet model (CSM) in hadroproduction 10 and the experimental data 
amounts to more than an order of magnitude and cannot be attributed to those the¬ 
oretical uncertainties arising from the ambiguities on the choice of a particular parton 
distribution function, the heavy quark mass or different energy scales. 

The colour-singlet model assumes that the heavy quark pair is produced in a colour- 
singlet state with the right quantum numbers already during the hard collision. The 
production cross section involves a perturbative calculation where the proper spin and 
angular momentum states are projected out, and the non-perturbative component of the 
hadronization comes from a single parameter to be identified with the wave function of 
the resonance at the origin (or its derivatives for non S'-wave states). The fact that the 
latter can be derived from potential models fitting the charmonia spectrum or from their 
leptonic decay modes, gives its predictive power to the colour-singlet model. However, 
as mentioned above, such predictions stand quite below the experimental yield at high- 
Pt, so that some modifications have to be done. In addition, it should be emphasized 
that the model is clearly incomplete since relativistic corrections due to the relative 
velocity of quarks in the bound state are ignored - expected to be sizeable especially 
for charmonia. Moreover, another hint of the incompleteness of the CSM comes from 
the existence of infrared divergences in the production cross section for P -wave states, 
which cannot be absorbed in its non-perturbative factor. 


In order to resolve all these difficulties, it has been recently argued that the heavy 
quark pair not necessarily has to be produced in a colour-singlet state at the short- 
distance partonic process itself |ll|. Alternatively, it can be produced in a colour-octet 
state evolving non-perturbatively into quarkonium in a specific final state with the 
correct quantum numbers according to some computable probabilities governed by the 
internal velocity of the heavy quark. This mechanism, usually named as the colour-octet 
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model (COM) []. can be cast into the rigorous framework of an effective non-relativistic 
theory for the strong interactions (NRQCD) deriving from first principles [|Hfj . 

However, the weakness of the COM lies in the fact that the non-perturbative param¬ 
eters characterizing the long-distance hadronization process beyond the colour-singlet 
contribution (i.e. the colour-octet matrix elements) are so far almost free parameters 
to be adjusted from the fit to experimental data, though expected to be mutually con¬ 
sistent according to the NRQCD power scaling rules, or to resort to non-perturbative 
models [|p| . 


On the other hand, an attractive feature of the colour-octet hypothesis consists of 
the universality of the NRQCD matrix elements (ME) entering in other charmonium 
production processes like photoproduction |jlbf. Let us look below in some detail at 


the way hadronization is folded with the partonic description of hadrons, focusing for 
concreteness on the couple of related papers 


16, 17 


In the first paper [JTB[ , Clio and Leibovich considered for charmonium production (in 
addition to the CSM) the 1 2 3 Si coloured intermediate state as a first approach, computing 
the matrix elements as products of perturbative parts for the short-distance partonic 
processes, and the colour-octet ME concerning the long-distance hadronization. Finally, 
a convolution of concrete parton distribution functions and the differential cross-section 
for the QQ production subprocess was performed, whereby the p t dependence of the 
charmonia production exclusively coming from the latter. 


In Ref. JTj 7 ] the same authors take into account further contributions from new 
coloured states (for more details see the quoted references) concluding finally that at 
high enough p t a two-parameter fit is actually required to explain the observed inclusive 
p t distribution of charmonia production at the Fermilab Tevatron. 


Notice, however, that such calculations were carried out based on a possibly oversim¬ 
plified picture of the hadronic interaction. Indeed, it is well-known for a long time that 
higher-order effects (K factors) play an important role in inclusive hadroproduction. 
In particular, beyond the primordial transverse momentum kt of partons in hadrons 
related to Fermi motion relevant at small p t , initial-state radiation of gluons by the 
interacting partons add up to yield an effective intrinsic transverse momentum which 

H 


certainly has to be considered in hadroproduction at high p t 


As we shall see, if 


overlooked at all, the effect on the fit parameters (and ultimately on the colour-octet 
ME’s) amounts to a systematic overestimate. 


In fact, one should distinguish the primordial transverse momentum of partons ow¬ 
ing to their Fermi motion (a non-perturbative effect) from the perturbative contribution 
dynamically generated via gluon radiation generally implemented in the event genera¬ 


tors by means of a parton shower algorithm 00- By effective k t we merge both effects 


1 For the sake of brevity in our theoretical introduction we do not expose another approach to heavy quarko- 
nia production in the same spirit as the COM but based on duality arguments known as the colour-evaporation 
model, referring the reader to j^] and references therein 

2 However, it is not yet completely clear how to separate uniquely both components from real experimental 

data [Q 
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under a common name, though the former is actually overshadowed by the latter at 
high p t . Of course, the PYTHIA treatment of the effective k t is not guaranteed to be 
perfect but, nevertheless, should give a reasonable estimate of such effects. 


In this work we have extended our earlier analysis of J/-0 hadroproduction at col- 
ZT| to the if}' state, affording a more complete test upon the validity of the 


liders 

NRQCD power scaling rules relating different long-distance matrix elements. The pa¬ 
per is organized as follows. In Section 2 we describe the hard processes considered for 
charmonia hadroproduction and we discuss their implementation in the event generator 
PYTHIA, described in more detail in the Appendix at the end of the paper where we 
also present a code for a fast simulation of charmonia using PYTHIA. In Section 3 we 
perform new fits to Tevatron data on prompt J/p and tj}' direct production extracting 
NRQCD matrix elements from them. We finally make some predictions for charmonia 
direct production at the LHC in Section 4. 


2 Implementation of the COM in PYTHIA 

Originally, the event generator PYTHIA 0 produces direct J/p’s and higher Xcj res¬ 
onances in hadron-hadron collisions according to the CSM. However, as already men¬ 
tioned this model fails off to account for the charmonia production rate at the Fermilab 
Tevatron by more than an order of magnitude. Consequently, we have implemented a 
code in PYTHIA in order to include the colour-octet mechanism for J /ip hadroproduc¬ 
tion ('0' generation has been based upon the same set of production channels) via the 
following cc 3 partonic processes: 

• 9 + 9 -> J/^ + 9 

• 9 + q -> J /P + q 

• q + q ^ J/i!> + g 


We have included in the simulation the 3 S[ 8 \ 1 S , q 8) and 3 Pj 8) (J = 0,1,2) states 
3 Si 


n 


m “£'i charmonia hadroproduction according to the COM 
tudes, taken from Refs. [00 are reproduced below. Firstly let us consider the 
intermediate state: 


The squared ampli- 

~s( 8) 


J 2 \ A ( 99 ^ J /^ 9)\ 2 = 

x 

x 


87r 3 a 3 27 (st + tu + us) — 19M 4 

~ 9M 3 [(s - M 2 )(t - M 2 )(u - M 2 )] 2 

[{t 2 T tu T u 2 ) 2 — M 2 (t + u)(2t 2 + tu + 2 u 2 ) + M 4 (P + tu + u 2 )\ 
<0|0 8 j/ ^( 3 A 1 )|0> (1) 


^2\A(gq^J/^q)\ 2 


167T 3 0! 3 
27 M 3 


s 2 + u 2 + 2 M 2 t 
su(t — M 2 ) 2 


[4(s 2 + u 2 ) - tu] < 0|0 8 W ( 3 A 1 )|0 > 


( 2 ) 
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128vr 3 ct 3 t 2 + u 2 + 2M 2 3 
81M 3 tu(3 -M 2 ) 2 


[A(i 2 + ii 2 ) - tu} < o|o 8 j/v, ( 3 5i) |o > 


( 3 ) 


where the barred summation symbol refers to average over initial and final spins and 
colours; 3, t and u stand for the Mandelstam variables of the short-distance subpro¬ 
cesses. We shall set M = 2 M c for both J/if and if' resonances in accordance with the 
assumption on the non-relativistic nature of heavy quarkonium made in the theoretical 
calculation of these amplitudes |16], [17|], i.e. vanishing relative momentum of the quarks 


in the bound state. 

The corresponding expressions for the 1 Sf' ) contributions are: 


J2\ A (99^J/^g)\ 2 = 

X 


Y^A(gq-^ J /il>q) | 2 = 

^2\A{qq^ J /i/jq)\ 2 = 


207r 3 a 3 ( 
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With regard to the 3 Pg contributions, they display altogether a similar 

(8) 


i.e. 


de- 

Hl- 


generate) transverse momentum behaviour as the 1 S'q component for p t > 5 GeV 
Thus from a pragmatic point of view, the generation of charnronia via intermediate 
P -wave coloured states becomes superfluous although they must properly be taken into 
account in the computation of the overall cross section as explained in more detail in 
the Appendix. 


We find from our simulation that gluon-gluon scattering actually stands for the dom¬ 
inant process as expected, the gluon-quark scattering contributes appreciably ~ 20%, 
whereas the quark-antiquark scattering represents a tiny contribution. Therefore we 
shall omit hereafter any further reference to the qq production channel. 


Using the same numerical values for the colour-octet matrix elements as reported 
in Tables 1 and II of Ref. [O] , if initial- and final-state radiation are turned off, there 


is a good agreement between the theoretical curve and the experimental points for 
both J/if and if>' (see Fig. 1), as should reasonably be expected . However, if initial- 
and final-state radiation [] are switched on, the predicted curve stands well above the 
experimental data, in accordance with the expected /q-kick caused by the effective 


,! It should be noted that initial-state radiation and final-state radiation have opposite effects in the pt spec¬ 
trum, the former enhancing the high pt tail whereas the latter softens the distribution. Indeed, in considering 
the process gg^J/ ip g in PYTHIA only the gluon evolves in the final state, though the energy (and momentum) 
of charmonium is modified as a consequence of the final-state machinery |i| 
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Figure 1: Curves obtained from PYTHIA (not fit) including the colour-octet mechanism 
for prompt J/'i) (left) and if' (right) production at the Tevatron using the same parameters 
as in Ref. [17]. The solid line corresponds to initial- and final-state radiation off and the 
dotted line to initial- and final-state radiation on. The MRSDO parton distribution function 
and M c = 1.48 GeV were employed as in [17]. All curves are multiplied by the corresponding 
muon branching fraction of charmonium. 


intrinsic transverse momentum of partons [18]. Accordingly, keeping radiation effects 


in the theoretical analysis it turns out that the values for the colour-octet ME’s have 
to be lowered by significant factors. 


For the sake of clarity it is important to point out that the inclusion in the event 
generation of multi-gluon emission by itself almost does not change the overall cross- 
section for charmonia production obtained from PYTHIA. (There is however an indirect 
bias effect slightly increasing the cross section cr(\rjj/^\ < 0.6) with radiation on relative 
to radiation off because of the kinematical cut on the charmonia pseudorapidity.) The 
main effect remains in the smearing of the transverse momentum distribution leading 
to the enhancement of the liigh-pi tail. 


3 Extraction of NRQCD matrix elements from Tevatron data 


In order to assess the importance of the effective intrinsic transverse momentum of 
partons we have made three different choices for the proton PDF Q: 


a) the leading order CTEQ 2L (by default in PYTHIA 5.7) 


b) the next to leading order MRSDO (the same as used in ||T7|] ) 


c) the next to leading order GRV 94 HO 

4 See [Q for technical details about the package of Parton Density Functions available at the CERN Program 
Library. References therein 
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B(J/V—>mV ) da(ppH> J/^+X)/dp, (nb/GeV) 


As mentioned before, the theoretical curves for the inclusive p t distribution of prompt 
J/ip and tp' stand in all cases above Tevatron experimental points if the set of parameters 
from |TjJ are blindly employed in the PYTHIA generation with initial-state radiation 
on. Motivated by this systematic discrepancy, we have performed new fits for prompt 
J/ip and ip' direct production at the Tevatron (feed-down to J/ip from radiative decay 
of Xcj resonances was experimentally removed). 

Those colour-octet matrix elements extracted from our analysis on J/ip production 
are reported in Table 1. In Fig. 2a we plot the resulting fit to Tevatron data for the 
CTEQ PDF, showing separately the three different components taken into account in 
this work. As a figure of merit for the overall fit we found x 2 /Ndf = 1-8. In the case of 
ip' production, the respective ME’s are reported in Table 2 and the corresponding plot 
for the CTEQ PDF is shown in figure 2b, with x 2 /Ndf = 0.4. 




Figure 2: Two-parameter fits to the experimental Tevatron data using CSM + COM, where 
initial- and final-state radiation were incorporated via PYTHIA generation, a) (left) J/ip 
and b) (right) ip 1 production. The different contributions are shown separately i) dotted line: 
CSM, ii) dashed line: 1 S'J 8) + 3 pj 8) , iii) dot-dashed line: iv) solid line: all contributions. 


Table 1: Colour-octet matrix elements (in units of GeV 3 ) from the best y 2 fit to Tevatron 
data on prompt J/ip production for different parton distribution functions. The error bars 
are statistical only. For comparison we quote the values given in Ref. [17]: (6.6T2.1)x 10 -3 
and (2.2T0.5) x 10 -2 respectively. 


matrix element: 

< 0|O^( 3 S' 1 )|0 > 

<0|0« /,/ '( 3 Po)|0> <o|o^ /I/, ( 1 s 0 )|o> 

M-? 1 3 

CTEQ2L 

(3.3T0.5) xlO -3 

(4.8T0.7) xlO -3 

MRSDO 

(2.1±0.5)xl0~ 3 

(4.4T0.7) xl0~ 3 

GRV 94 HO 

(3.4T0.4) xlO -3 

(2.0±0.4) xlO -3 
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Table 2: Colour-octet matrix elements (in units of GeV 3 ) from the best y 2 fit to Tevatron 
data on prompt ip' production for different parton distribution functions. The error bars are 
statistical only. For comparison we quote the values given in Ref. [17]: (4.6T1.0)x 10 -3 and 
(5.9±1.9)xl0~ 3 respectively. 


matrix element: 

< o|of psoio > 

<o|of ( 3 P o )|0> , <o|of ( 1 s 0 )|o> 

M2 1 3 

CTEQ2L 

(1.4±0.3)xl0 -3 

(l.l±0.3)xl0 -3 

MRSDO 

(l.l±0.3)xl0~ 3 

(0.94±0.23)xl0~ 3 

GRV 94 HO 

(1.3±0.2)xl0~ 3 

(0.12±0.15)xl0~ 3 


Notice that all the ME’s are lowered systematically with respect to the values re¬ 
ported in Ref. [[17] for both J/ip and ip' []. Basically this is because initial-state radiation 
enhances the differential cross section for p t > 4 GeV, so that the non-perturbative pa¬ 
rameters have to be re-adjusted to restore the fit to the experimental points. Let us also 
remark that the combined 1 Sq'\ 3 P 0 ® ME’s are lowered by larger factors than those for 
the 3 S® contribution. Consequently, the power counting rules derived from NRQCD 

0 


are more faithfully satisfied in the overall than in Ref. 


all of them should scale 


like Mpv'. It is especially significant that the numerical values obtained in our analysis 
for the 1 4 8) + 3 P 0 (8) ME’s are smaller by an order of magnitude with respect to those 
found in Ref. |17| in accordance with the common belief widely expressed in literature 
q [24|| . In this sense, theoretical expectations from the COM might be reconciled with 


HERA data on J/ip photoproduction [|H| 


The contributions from the colour-octet 3 S® and 3 S® + 3 P® states to the differ¬ 
ential cross section for J/ip and ip’ production are shown separately in Fig. 2. Notice 
that the 1 S^ an d 3 P® components dominate in both cases over the other production 
channels at moderate pt, whereas the 3 S'® channel takes over at large pt becoming 
asymptotically responsible of the full production rate. 

Let us remark the especially relative smaller numerical values of the long-distance 
1 S'®, 3 P® matrix elements regarding the GRV PDF. This is in a qualitative agreement 
with the fact that the MRSDO and CTEQ2L parton distribution functions give smaller 
values for the gluon distribution function at small partonic x than the GRV PDF (the 
gluon density of the latter is steeper at small x ). Hence, a gluon density with a steeper 
behaviour at small x steepens the p t - spectrum for the 3 S'® component amounting to 
a larger production rate at moderate p f , decreasing the 1 S'® + 3 P® contribution from 
the combined fit to the experimental inclusive distribution. 

5 A caveat is in order however. In our generation we have not included the leading log QCD corrections into 
the colour-octet cross section (causing the depletion of the frag mentation function of gluons into charmonium at 
high-pt) through the interpolating formula (2.22) of Ref. [ha. However, final state radiation was incorporated 
in our analysis in contrast to Ref. Bill: softening the pt inclusive distribution as mentioned in footnote #3. 
The neglect of those higher order QCD effects should mainly affect the high-pt region, hence the 3 Si matrix 
element which could be somewhat underestimated. 
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4 Charmonia Production at the LHC 

Reliable calculations of inclusive production rates of charmonia are of considerable inter¬ 
est at LHC experiments in many respects 0,0, as briefly outlined in the Intorduction. 
Therefore, an order-of-magnitude estimate of its production rate becomes very valuable 
at the present stage of the LHC experiments. 




Figure 3: (left) Different contributions to J/if hadroproduction using the CTEQ PDF. The 
curves are labelled in the same way as in Fig. 2. ; (right) Our prediction for prompt J/if 
direct production at the LHC using PYTHIA with the colour-octet matrix elements from 
Table 1: a) dotted line: CTEQ 2L, b) dot-dashed line: MRSDO, and c) solid line: GRV HO. 
A rapidity cut \y\ < 2.5 was required. 


Keeping in mind such a variety of interests, we have generated direct J/if 's and if's 
in proton-proton collisions at LHC energies (center-of-mass energy = 14 TeV) by means 
of our “modified” version of PYTHIA with the colour-octet ME’s as shown in Tables 
1 and 2 - i.e. after normalization to Tevatron data. The different contributions to J/if 
hadroproduction are depicted altogether in Fig. 3 for the CTEQ PDF. 


Figure 3 also illustrates the p t inclusive distributions of direct prompt J/if's at 
the LHC obtained for the three PDF’s employed in our study. Comparing them with 
the distribution obtained by Sridhar [^4|, we find our predictions standing below his 


theoretical curve at large pt . In fact this is not surprising since Sridhar considered only 
the 3 S[ S) coloured intermediate state, whose ME was assumed larger than the one used 
in this work. The argument ends by noticing that at high enough p t the dominant 


production comes from the 3 S[ 8 ' 1 
falls off faster than the former. 


channel as the combined 1 Sn 8 ^ + 3 Pl 8) contribution 


In Fig. 4 we plot our prediction for direct if' production at the LHC using the 
CTEQ PDF, together with the J/if curve. From the figure it is clear that muon pair 
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Figure 4: J/ijj (solid line) and ip' (dash-dotted line) direct production at the LHC according 
to the colour-octet ME’s reported in Tables 1 and 2 using the CTEQ PDF. The cut \y\ < 2.5 
was required in both cases. Both curves are multiplied by the respective branching fraction 
of charmonium into muons. 


production from direct J/ip is larger by more than an order of magnitude from J/ip as 
compared to ip 1 . It should be mentioned that our prediction for the ip' production rate 
differs substancially with respect to the corresponding theoretical curve by Sridhar |24 


all along the p t range - by an order of magnitude even at moderate p t in contrast to the 
J/ip case. 


5 Conclusions 


In this paper we report on the extraction of long-distance colour-octet matrix elements 
in a Monte Carlo framework from charmonia production at hadron colliders, already 
studied in our previous works [^T|, extending our analysis to the ip' resonance. We 
explain in more detail the procedure followed in our investigation, also providing a code 
for a “fast” simulation of charmonium in hadron collisions based on the colour-octet 
model and compatible with the Tevatron results. 


All three possible lowest order short-distance processes contributing to J/'ijj hadropro- 
duction have been taken into account, namely, g + g^-J/ip + g, g + q—>J/ip + q 
and q + q^J/ip + g, contributing at order af to charmonia production including 
the colour-octet mechanism. The same set of channels have been considered for ip' 
hadropro duct ion. 


We have investigated higher-order effects induced by initial-state radiation, con¬ 
cluding that the overlooking of initial radiation leads systematically to a significant 
overestimate of the fit parameters. The new colour-octet matrix elements obtained for 
three different PDF’s from Tevatron data are shown in Tables 1 and 2, which have to be 
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regarded as reasonable estimates due to the experimental and theoretical uncertainties. 
In particular, the numerical values for the combined 1 5 'q 8 ' ) + 3 Pq 8 ' ) matrix elements are 

and might resolve 

H- 


smaller by an order of magnitude than those obtained in Ref. [17 


the difficulties of the COM found in J/-0 photoproduction at HERA 
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Appendix 

A Code for Implementing the Colour Mechanism in PYTHIA 


We have modified the PYTHIA routine PYSIGH (ISUB=86) || for the partonic pro¬ 
cess: 


g + g -»■ JH> + g 


including those contributions from the 3 S[ 8 \ 1 S'q 8) anu i j 
to the expressions (1) and (4). The ME’s of the latter components are related through 
the heavy quark spin symmetry at leading order in 1 /M c : m 


1 S'n 8) and 3 P ( , 8) coloured states according 


< O s ( 3 Pj ) > = (2 J + 1) < 0 8 ( 3 Pq) > 


Moreover, it should be stressed that the colour-octet ME’s associated to the 1 S and 
the 3 Pj 8 ^ channels cannot be disentangled from Tevatron experimental data since both 
are degenerate in the sense that they display a similar p t behaviour for pt > 5 GeV 
17| . All three P-wave contributions altogether amount to about three times the 1 S'o 
contribution when setting < O 8 ( 3 P 0 ) >= M% < O 8 ( 1 S' 0 ) >. Therefore one may write 
the following approximate equality: 


So) > < Q 8 ( 3 Pq) > 
Ml 
(AT) 

and only the 1 S^ component has to be simulated in practice. 


da 

dpt 


[‘sn + E 4r[ 3 d 8) ] - 3 x 


da 


J= 0,1,2 


dp t 


< O 8 ( 1 S 0 ) > dp t 


'SP} 


x 


< Os 



Figure 5: Histograms for ^~)xda/dpt{pp^J/'tjjX) obtained from PYTHIA 

using the CTEQ PDF at 14 TeV center-of-mass energy. Solid line: all contributions; dashed 
line: gg contribution. 
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On the other hand, we have generated with PYTHIA an event sample via the gluon- 
quark subprocess g + q-^Jftp + q according to the expressions (2) and (5), “adapting” 
for this end the PYSIGH routine ISUB=112. We found that this production channel 
amounts to a fraction of about 20% of the total yield of direct J ftp's, almost constant 
over the p t range under examination (see Fig. 5). Consequently, one may reduce 
considerably CPU time in the event generation by a slight change on the gg-^>J/tpg 
routine, including gq scattering by means of a scaling factor. Let us note however 
that all the fits and plots shown in this work were performed using the gq event sample 
obtained from PYTHIA as explained above. The routine presented below permits a fast 
generation of direct Jftp's (or tp 1 with some few changes in the values of the parameters 
and variables such as PARP(38), SQM3, ...) in hadronic collisions at the LHC. 

Thus, we have rewritten the ISUB=86 routine of PYTHIA as: 


ELSEIF(ISUB.EQ.86) THEN 
C...g + g -> J/Psi + g. 

C 

C... Color-Octet Matrix Elements for J/psi 
C (CTEQ 2L) 

C 

SQM3=(2.*1.48)**2 

X083S1=0.0033 

X081S3P=0.0048 
C 

C COLOR SINGLET 
C 

FACQQG1=C0MFAC*AS**3*(5./9.)*PARP(38)*SQRT(SQM3)* 

& (((SH*(SH-SQM3))**2+(TH*(TH-SQM3))**2+(UH*(UH-SQM3))**2)/ 

& ((TH-SQM3)*(UH-SQM3))**2)/(SH-SQM3)**2 

C 

C COLOR OCTET (3-S-l) 

C 

FACQQG2=-C0MFAC*AS**3*(1./18.)*3.1416/SQRT(SQM3)* 

& (27.*(SH*TH+TH*UH+UH*SH)-19.*SQM3**2)*((TH**2+TH*UH+UH**2)**2- 

& SQM3*(TH+UH)*(2.*TH**2+TH*UH+2.*UH**2)+SQM3**2* 

& (TH**2+TH*UH+UH**2))*X083S1/ 

& ((TH-SQM3)*(UH-SQM3))**2/(SH-SQM3)**2 

C 

C COLOR OCTET (l-S-0)/3 + (3-P-O) !!! 

C 

FACQQG3=1000.*C0MFAC*(5.*3.1416*AS**3/(4.*SQRT(SQM3)))/SH* 

& (SH**2*(SH-SQM3)**2+SH*TH*UH*(SQM3-2.*SH)+(TH*UH)**2)/(TH*UH)* 

& ((SH**2-SQM3*SH+SQM3**2)**2- 

& TH*UH*(2.*TH**2+3.*TH*UH+2.*UH**2))*X081S3P/ 
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& ((TH-SQM3)*(UH-SQM3))**2/(SH-SQM3)**2 

FACQQG3=FACQQG3/1000. 

C 

C ... Scaling factor for g + q -> J/psi + q 
C 

SCALGQ=1.25 
C 

IF(KFAC(1,21)*KFAC(2,21).NE.O) THEN 
NCHN=NCHN+1 
ISIG(NCHN,1)=21 
ISIG(NCHN,2)=21 
ISIG(NCHN,3)=1 
SIGH(NCHN)=FACQQG1 
NCHN=NCHN+1 
ISIG(NCHN,1)=21 
ISIG(NCHN,2)=21 
ISIG(NCHN,3)=1 
SIGH(NCHN)=FACQQG2*SCALGQ 
NCHN=NCHN+1 
ISIG(NCHN,1)=21 
ISIG(NCHN,2)=21 
ISIG(NCHN,3)=1 

SIGH(NCHN)=3.*FACQQG3*SCALGQ 

END IF 


The following variables, originally not present in PYTHIA, have been specifically 
introduced in our routine: 

• X083S1: colour-octet ME associated to the 3 S'i component 

• X081S3P: linear combination < O 8 ( 1 S , 0 ) > /3+< O 8 ( 3 P 0 ) > /Af? 

• SCALGQ: scaling factor taking into account gq production 

A lower cut-off was used throughout in the event generation (ptmin = 1 GeV by de¬ 
fault in PYTHIA for processes singular at p*—>0) thereby avoiding the problematic low- 
Pt region for charmonia hadroproduction. The effect of such a cut on the experimentally 
observed pt range ( p t > 4 GeV) is rather small, representing a slight underestimation 
of the production rate in the event generation with radiation on as a consequence of 
the ISR effect [Pf] . 

Lastly, let us point out that the theoretical curves appearing in all our plots were 
obtained from histograms (filled with the corresponding PYTHIA ouput data) fitted 
by means of the four-parameter p t function: 

Pt 2 

F[oti, ct 2 , a 3 , aq; p t ] = aq --* (A.2) 

(«3 +Pt) 4 
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providing an excellent fit in the whole p t range examined. 
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